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This article presents some results in terms of the performance improvement of 
a multi- feed array configuration over the usual single feed system when an adaptive 
least-squares algorithm is applied for the signal reconstruction . The article presents 
two novel versions of the least-squares algorithm , one of which is based on the max- 
imization of the signal-to-noise ratio while the other is based on the deconvolution 
of the received signal field. These algorithms have been developed for the purpose 
of minimizing degradations arising from various sources, which can severely limit 
the performance (gain) of a single- feed system. 


I. Introduction 

The development of multi-element array process- 
ing techniques has many potential applications for the 
Deep Space Network (DSN). These include signal recon- 
struction for both X-band and Ka-band communications 
[1-4] and electronic pointing to augment existing mechan- 
ical pointing techniques. For all of these applications, 
multi-element array processing can generally provide sig- 
nificant performance improvements over single-feed an- 
tenna configurations, which are predominantly used in the 
DSN. 

This article presents some results in terms of the per- 
formance improvement provided by a linear multi-feed ar- 
ray incorporating the proposed adaptive least-squares al- 
gorithm over a single-feed array system, as applied to the 
signal reconstruction problem. While the assumed linear 
array geometry is idealized, the results of this analysis pro- 


vide an indication of performance improvements that can 
be achieved with adaptive, multi-element array processing. 
This article describes two versions of the least-squares al- 
gorithm, one of which is based on the maximization of the 
signal-to-noise ratio while the other is based on deconvo- 
lution of the received signal field. Here, instead of trying 
to model the signal degradations in terms of deterministic 
equations in evaluating the performance of the algorithm, 
it is assumed that these degradations are “unknown” to 
the algorithm and vary with time. The algorithm tries to 
implicitly estimate these degradations in an adaptive man- 
ner from the samples of the noisy received signal. On the 
basis of these measurements, it computes a set of weights 
for combining signals at the outputs of various feeds in or- 
der to maximize the signal-to-noise ratio of the combined 
signal. 

For the purposes of illustrating the basic concepts in- 
volved with adaptive array processing, this article presents 
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the results for a 16-element linear-array feed system. The 
performance of the least-squares algorithm is to a first or- 
der determined by the signal- to- noise ratio of the received 
signal, the number of feeds in the configuration, and the 
time constants at which the received signal field is varying 
with time. 

In practice, it may be possible to simultaneously cor- 
rect for multiple degradations arising from different sources 
and having different time constants. These degradations 
may be induced by wind, gravitational loading, or antenna 
pointing errors. Simultaneous correction for such degrada- 
tions could be achieved by adjusting the time constants of 
the algorithm to track the fastest mode, in which case the 
slower modes would be estimated sub-optimally. Alterna- 
tively, one could track the most significant mode, thereby 
essentially ignoring the faster but less significant modes. 
More sophisticated techniques could also be used for sep- 
aration of the modes and tracking of them separately. 

II' Array Configuration 

The specific array configuration of interest in this ar- 
ticle corresponds to a linear-feed array distributed across 
the focal plane of an antenna. The array outputs are then 
fed to a parallel receiver bank as indicated in Fig. 1. As 
shown in the schematic diagram of Fig. 1, the signal out- 
puts from the feed elements are amplified by r f amplifiers. 
Assuming that all of the amplifiers have equal gain and 
noise temperature, the output of the zth amplifier can be 
written as 

rj(<) = ylj(i)cos(w c < + + rii(t ) (1) 

where Aj(2) and @i(t) are the signal amplitudes and phases, 
u c denotes the signal carrier frequency, and n*(t) is a zero- 
mean white Gaussian noise of one-sided spectral density 
Nq. The noise is also assumed to be spatially uncorrelated, 
i.e., E [nj(2)nj(t)] — 0 for i / j and i,j = 1,2,.., N, Under 
ideal conditions and assuming a plane-wave normally inci- 
dent on the antenna aperture, the amplitude of the center 
feed would be equal to y/2 P ( P denotes the normalized 
power received by the antenna), while the remaining feeds 
will have nearly zero amplitude. However, array degrada- 
tions can disperse the signal amplitude (and phase) spa- 
tially over N feeds. These degradations can arise due to 
various sources, such as gravity, thermal fields, wind, and 
atmospheric turbulence. 

In the presence of such degradations, the adaptive sig- 
nal processing algorithm then combines the N feed outputs 
in a coherent manner so as to optimise some performance 


index, such as the signal- to-noise power ratio of the com- 
bined signal. This processing can occur either at r f or 
can be equivalently done at the baseband. Alternatively, 
in some possible imaging applications, it may be desired 
to reconstruct the complete focal plane field, i.e., obtain 
N output signals that are close to the outputs of the focal 
plane fields in the ideal antenna case. 1 Note that in the 
more general case, more than one feed may have signifi- 
cant amplitude if the source is not a point source and the 
antenna is capable of resolving such a composite source. 
In the limiting case, one may simply use the center out- 
put of the reconstructed field and ignore the others, thus 
achieving an alternative combination of the input signals. 

For the purposes of signal processing, the N received 
r f signals r*(t); i = 1,2, are down-converted and 

quadrature sampled to obtain the sampled version of the 
complex baseband envelope gi(t) of the r f signal r*(t) with 

r i(t) = Re{g i (i)e^ t } 

9i (t) = Mtyw + Vi (t) ( 2 ) 

9i(t) = {r<(0 + m}e~ jUct 

In Eq. (2) above, r*(i) denotes the Hilbert transform 
of r*(i), and V{{t) is the complex envelope of the bandpass 
noise n,(£). 

III. Signal Combining Via Adaptive 
Least-Squares Algorithm I 

As shown in Fig. 1, the adaptive algorithm deter- 
mines the time-varying complex-valued weighting coeffi- 
cients w \ (&), ...,^tv(^) on the basis of signal samples 
i — 1,2, ..., N\ and j = 1, 2, ..., k according to some appro- 
priate optimization criterion. The algorithm is adaptive 
in the sense that if the signal amplitudes and phases (A,- 
and Qi) remain relatively constant with time, then with 
increasing value of k , the algorithm achieves increasingly 
accurate estimates of these parameters, and the weighting 
coefficients con% r erge asymptotically to their theoretically 
optimum values with an exponential convergence rate. On 
the other hand, if these parameters are time-varying, then 


1 V. Vilnrotter, “Ka-Band Array Signal Processing Progress 
Report,” JPL Interoffice Memorandum No. 331-88.5-047 (in- 
ternal document), Jet Propulsion Laboratory, Pasadena, Cal- 
ifornia, November 1988. 
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the algorithm tracks these variations and the weighting co- 
efficients are truly time-varying (there is no tendency for 
Wi(k) to converge to some constant value). 

Denoting by w(k) and g(k) the weighting coefficient 
vector [wi(k)w 2 (k) . . .WN(k)] f and the measurement vec- 
tor [gi(k)g 2 (k) . . . <7jv(&)] / respectively, then the familiar 
least-square optimization criterion is to select w(k) so as 
to minimize the following index 


One may note that the entity to be inverted in Eq. (7) is 
only a scalar. Decomposing the sum in Eq. (5) as 

+ 9_(k) 

i = i 

and substituting Eq. (7) for P(fc), we obtain the following 
expression for w.ELs(k)- 


Jk = y:i i -^(jfe)£(i)i 2 

3 = 1 


( 3 ) 


&ELs( k ) — £(* — 1) 



with respect to the weight vector w(k) for k = 1,2, .... In 
the above, the superscript H denotes conjugate transpose 
while / represents just the transpose of a matrix. The op- 
timal solution, termed least-squares estimate of w(£), is 
given by (assuming k > N) 


WLs( k ) = 

If the distortion process is time-varying, then it is more 
appropriate to replace the index J* by the one obtained 
by multiplying the summand in Eq. (3) by for some 
0 < A < 1, minimization of which yields the following ex- 
ponentially data- weighted least-squares estimate for w(k). 

MLELs(k) = ji; X k - j 9_(j)9_ H (j) j E A (5) 

Note that in the adaptive algorithm’s present non-recursive 
form, Eq. (4), it is required to invert an (N x N) matrix for 
every time instance k in the computation of w(k), which 
is somewhat computationally intensive. This problem can 
be overcome by replacing the estimate in Eq. (5) with its 
recursive form, which is obtained as follows. 

Denoting by P(fc) the matrix inverse in Eq. (5), then 
the matrix PT l (k) has the following update. 


Ei0V*0')f E &O0 ( 4 ) 

j = l J J = 1 


- [A + £ "(fc)P(fc-l)£(fc)]" 
x £(Jfc - l)g(k)g H (k)P(k - 1) 

x |E Ai_ 1 _i £ 0 )| + ™#) ( 8 ) 

By noting that the first term in Eq. (8) and the product 
of the last two factors in the second term both are equal 
to UlELs(k “ 1)» Eq. (8) may be rewritten as 

mELs(k) ~ HlELsi^ ~ 1 ) + £(*)#) 

- [A + £ fl (fc)£(fc-l)£(fc)] _1 
x P(k - 1 )l(k)g n (k)w ELS (k - 1) (9) 

Post multiplying both sides of Eq. (7) by g(k) and with a 
simple algebraic manipulation, it follows that, 

[a +£ JI (*)£( i-l)£(*)]~ 1 £(*-l)£(*) =£(*)£(*) ( 10 ) 

With the substitution of Eq. (10) in the last term of 
Eq. (9), one obtains the recursive version of the algorithm 
given below. 


™.ELs( k ) = ™ELs(k - 1 ) + P( k )9_(k) 


£ -I (*) = A£ _1 (*-l) + £(*)£*(*); it = 1,2,... (6) 

Application of the matrix inversion lemma [5] to Eq. (6) 
yields the following desired recursion for P(fc). 

P{k) = \~ 1 {p(k-l)- [\ + g H { k )P{k - l)l(k)]~ X 

x £(Jk - !)£(*)£"(*)£(* - 1)} (7) 


x [l — 9_ H ( k )a±ELs( k ~ 1) 

£(*) = A" 1 (£(Jfc - 1) - [A + g H (k)P(k 

x P(k - 1 )g(k)g H (k)P(k - 1) j; 
* = 0 , 1 , 2 ,... 
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If the complex field g(k) is a wide-sense stationary pro- 
cess, a considerable simplification in computations may be 
achieved by replacing £(£) with an appropriate constant 
matrix in the first recursion of Eq. (11) and dropping the 
second recursion. 

IV. Maximization of Signal-to-Noise Ratio 
via Modified Least-Squares Algorithm I 

In some applications, it may be more appropriate to 
maximize the signal- to-noise power ratio at the combiner 
output. The noise variance at the output of the combiner is 
equal toJ^j |tD,| 2 £ , [|i/,‘| 2 ] = <r 2 ||2)|| 2 with a 2 denoting the 
variance of the sampled version of the complex baseband 
process i ^(t) in Eq. (2). Thus, as shown in the Appendix, 
an effective maximization of the output signal- to-noise ra- 
tio can be achieved by minimization of the index given in 
Eq. (3) subject to the equality 

INI 2 = K (12) 

for some constant A", or by simply minimizing the following 
index, 

Ell- t£ ff (^)£(i)| 2 + 7«(IN| 2 - K) (13) 

j = 1 

with respect to w and j(k ), where 7(fc) is the Lagrangian 
multiplier. Differentiation of Eq. (13) with respect to w 
yields the following constrained least-squares estimate for 
w in terms of j(k) as 

MctsW = |E£^)£ ff w + 7 ^) 7 | ( 14 ) 

Substituting Eq. (14) into Eq. (12) yields an equation for 
the unknown 7 (k) } which can be solved for 7(fc). Substi- 
tuting 7 (k) back in Eq. (14) provides the constrained op- 
timum solution for the weighting coefficient vector. Note, 
however, that there is no close-form solution for 7(fc) and, 
thus, some numerical optimization techniques may need 
to be applied to obtain w CLS . A simplified solution is 
obtained by selecting some appropriate value for 7 (Jfc) in 
Eq. (14) and then normalizing the estimate to have its 
norm equal to one. In an exponentially data- weighted ver- 
sion of Eq. (14), both the summands are multiplied by A* _J 
where A is the exponential data-weighting coefficient. With 
these modifications, Eq. (14) has the following equivalent 
form. 


Pr\k) = A P~\k - 1) + 70 / + g(k)l H (k) 

±(k) = \±(k - 1) + g{k) 

w(k) = £(*)£(*) (15) 

ULnls = 2£(*)/||w(fc)|| 

7 (*) = 7o(l-A*)/(l-A) 

Note that Eq. (15) requires a matrix inversion for each 
value of k for which w(k) is desired. An approximate re- 
cursive form for Eq. (15), which does not require matrix 
inversion, may also be derived by applying the matrix in- 
version lemma. Note that in Eq. (15), the higher value 
of 7o results in the higher relative weight assigned to the 
noise variance at the combiner output. The initial values 
for P _1 and at k = 0 may simply be selected equal to 
zero. 

V. Least Squares Algorithm II 

In an alternative solution, it may be assumed that 
the received focal plane signal is the result of the spatial 
convolution of the ideal signal (in the absence of any distor- 
tion) with an unknown filter response representing various 
distortions from all sources, including antenna surface de- 
formations due to gravity, wind, antenna pointing errors, 
turbulence, etc., i.e., 

g(k) = B_ 2£(&) 4* v(k) (16) 

where B is a Toeplitz matrix, AC is the focal plane sig- 
nal vector that under appropriate sampling is equal to 
[0 ... 010 ... 0] in the ideal case of plane- wave with no dis- 
tortion, and v_ is the additive noise vector. The matrix B_ 
includes any distortion effects, pointing errors, etc. For the 
linear array case under consideration, 5 can be approxi- 
mated by a circular matrix for large N or becomes identical 
to a circular matrix provided the vector X_(k) and the sig- 
nal vector g(k) are zero-padded. Thus, it is assumed that 
B is circular. In the absence of noise, it is observed that 
A" = g where P -1 is also a circular matrix. How- 
ever, B is unknown, and it needs to be estimated from the 
noisy observations. Or, more directly, B_~ l is estimated as 
follows. Rewriting the model Eq. (16) as 

X(k) = F*g(k) + t ;(*); i = 1, 2, . . . (17) 

where F is an unknown circular matrix to be estimated 
from the given measurement g(k ); k = 1,2 ... ,n. Letting 
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f T = [fi , / 2 . . . /n] denote the first row of the matrix 
and £^(ifc) denote the vector obtained by cyclically shifting 
g(k) right i times, the following equivalent signal model is 
obtained. 


’o' 


£ T (*) ' 


7r 


Vl(fc) 

0 









£[(*) 


n 


{k) 

1 

= 




+ 


0 


* 





0 




Jn, 


v^f(k) _ 


(18) 

A least-squares estimate for / can be obtained from the 
signal model Eq. (18) and in its non-recursive form is given 
by 

/(*) = jll £*(■*’)£?(•>)) j E£n/ 2 0') ( 19 ) 

A recursive form for the estimate f(k) may also be derived 
following the steps used in obtaining Eq. (11). Here re- 
cursion is over both the signal sample vector g(k) and its 
circular shifts. With appropriate initial values P(0,0) and 
estimate /(0,0), one has the following recursion in terms 
of the indices k and with k denoting time and j denoting 
the cyclic shift of the received signal vector, 

P{k,j) = x- 1 {p(k,j- 1) 

xP(k, j - 1 fy.xWg-tW&’i - !)}; 

3- 1,2,... 

P(k + 1,0) = P(k,N)-, k= 1,2,... 

/(*, j) = i{k,j-l) + P(k, (*) 

x -£^. ,(*)/(*> j -1 )]; 

3 = 1,2.... 

f(k + 1,0) = f(k,N); k= 1,2,... 

( 20 ) 


where 

& = 1, 3 = [N/2\ 

= 0, j±\Nf 2J 

where [x\ denotes the least integer greater than or equal 
to x for any real x. The circular spatial convolution of 
/* with the received signal g(k) yields the reconstructed 
signal vector /i(fc). In the case of perfect reconstruction 
(deconvolution), the central element of the vector h is the 
combined signal, while the remaining elements would be 
zero. 

VI. Simulations 

The performance of the least-squares algorithms of 
the previous section is presented here in terms of simu- 
lations. In the following simulations, a linear feed array 
is considered. For the purposes of these simulations, the 
received signal focal field is generated by spatially Fourier 
transforming a simulated linear aperture plane array. The 
signals in the simulated aperture are assumed to be of 
equal amplitude but with completely independent phase 
processes. Also, for the purposes of simulations, each of 
these phase processes is assumed to be a moving aver- 
age process of a specified correlation interval I\ and vari- 
ance (steady-state). The effectiveness of the different least- 
squares algorithms is measured in terms of the power ratio 
(in dB) of the reconstructed (combined) signal to the total 
received power (which would be concentrated in the cen- 
tral feed element under ideal conditions). In addition, the 
performance of the least-squares multi-element combining 
algorithm is compared against traditional single-element 
processing. It should be noted that for all of these sim- 
ulation experiments, the noise variance at the combiner 
output is matched to the noise variance at the output of 
a single feed (the weighting parameter vector w is normal- 
ized to have unit norm). 

In the simulations reported here, it is assumed that 
the received signal is unmodulated with the signal field 
amplitude equal to 1 in the aperture plane, corresponding 
to a signal amplitude Aq equal to 16 at the center feed in 
the focal plane (ideal case). In case of data modulation, a 
decision-directed approach may be used to remove the data 
modulation from the received signal. For the purposes of 
simulations, the sampling period is normalized to 1 sec. 
Thus, the variance cr 2 of the sampled complex envelope of 
the noise at any of the feed outputs used in the simulations 
is given by 

<r 2 = A 2 0 /N 0 f s (21) 
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where f» is the sampling frequency; No is the one-sided 
noise spectral density; P — A^/2 is the received-signal 
power; and the algorithm’s performance is plotted as a 
function of the carrier-to-noise spectral density ratio: CNR 
= (P/No) (in dB-Hz). 

In Fig. 2, the reconstructed signal amplitude and 
phase estimates as a function of time measured in num- 
ber of samples for the least-squares algorithm, for a CNR 
of 10 dB-Hz are plotted. It may be observed from Fig. 2(a) 
that there is a considerable signal loss compared to the 
received signal amplitude equal to 16 for the ideal case. 
In the simulations, the initial estimate for the parameter 
vector w is selected to be [a a . . . a] with a = 1/16. As 
is apparent from Fig. 2(b), the phase-estimation error of 
the reconstructed signal is much smaller compared to the 
rms phase error of 1 rad introduced in the simulated re- 
ceived signal field. Figure 3 plots the corresponding results 
for the modified least-squares algorithm for the same set 
of signal parameters and for To = 100. Comparison of 
Figs. 2(a) and 3(a) shows a very significant performance 
improvement due to the proposed modification of the least- 
squares algorithm. Results are plotted in Fig. 4 for the case 
of 20 dB CNR with the parameter To equal to 200. For this 
case, the signal amplitude and phase of the reconstructed 
signal are quite close to their respective values for the ideal 
case. The signal amplitude loss for this case is only 1.09 dB 
relative to the ideal case, and the rms phase error (after 
adaptive combining) is 0.1 rad. 

A. Simulation Results for Least-Squares 
Algorithm I 

Figure 5(a) plots the sample estimates of the signal 
power loss (compared to the ideal case) for the standard 
least-squares algorithm I. The signal loss Pi is simply com- 
puted as Pi = 20Iog 10 (A rm */16) with 


A r ms — 


1 


M 




where 0, is the phase of the combined signal at the ith 
sampling instance. It is apparent from these figures that 
although the least-squares algorithm is optimal with re- 
spect to the prediction error criterion, it is not satisfactory 
in terms of the signal-to-noise ratio of the combined sig- 
nal. For the case of simulated phase dynamics, there is an 
asymptotic signal loss of 4 dB for the high CNR (~ 20 dB) 
case. 


Figure 6 plots the performance of the modified least- 
squares algorithm with To — 100 and with the same set 
of signal parameters as for the case of Fig. 5. The results 
are computed for three different values of the weighting 
coefficient: A equal to 0.925, 0.95, and 0.975. Compari- 
son of Figs. 5(a) and 6 shows a dramatic improvement in 
performance due to the proposed modification. Thus, the 
asymptotic signal loss for this case and with A = 0.925 is 
only 1.3 dB compared to a 4 dB loss for the standard least- 
squares algorithm. Increasing the value of To or reducing 
A may further reduce the signal loss. 

It is not difficult to understand this marked perfor- 
mance difference between the two algorithms. By mini- 
mizing the sum of the norm square of the estimated sig- 
nal error and the noise variance of the combined signal 
output, the standard least-squares algorithm produces a 
relatively large noise variance at the output. This is so 
because at high CNR the contribution of noise to the to- 
tal error is relatively much smaller than the contribution 
due to signal error (measured as a fraction of total sig- 
nal), and, thus, the former is essentially ignored by the 
algorithm. In the constrained least-squares algorithm, by 
increasing the relative weighting attached to the noise vari- 
ance, the signal-to-noise ratio at the combiner output is 
increased. In fact, as is shown in the Appendix, whereas 
the standard least-squares algorithm effectively maximizes 
the signal plus noise power output (subject to a near- 
orthogonality constraint), the constrained algorithm ac- 
tually maximizes the signal-to-noise power ratio (subject 
to a similar near-orthogonality constraint) at the combiner 
output. 


where Ai is the reconstructed signal amplitude at the ?*th 
sampling instance, and M is the number of sample values 
selected to be equal to 200 for these simulations. Fig- 
ure 5(b) plots the signal estimation rms phase error 0 rm5 
computed as 


The signal loss results presented in Figs. 6(a)-(c) can 
also be used to infer the performance gains provided by 
the multi-element least-squares technique over center-feed 
processing. In particular, we can define the array process- 
ing gain as the difference (in dB) between the combining 
losses of the least-squares and center-feed curves. These 
data are plotted in Fig. 6(d). 


© 


rm« 



M 


E »? 


As seen in Fig. 6(d), the array processing gain is a 
function of both CNR and A. Best results occur for large 
CNR (> 15 dB) and low values of A (A = 0.925). In par- 
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ticular, a maximum array gain of 3 dB (Fig. 6(a)) can be 
achieved for this simulation example. For lower values of 
CNR (below 0 dB), the estimation variance inherent in 
the least-squares algorithm degrades the array processing 
gain to such an extent that center-feed processing actually 
provides better performance. 

However, it must be stressed that the simulation 
example considered here corresponds to a relatively 
high dynamic scenario, e.g., as compared with typical 
mechanically-induced array degradations. Lower dynamic 
scenarios permit longer time constants for the least-squares 
algorithm (A — ► 1 and/or longer sampling period), thereby 
reducing the algorithm estimation variance. Thus, for low 
dynamic scenarios, positive array gains can be achieved 
over a wider range of CNRs (including CNR <C 0 dB) by 
utilizing time constants that are essentially matched to 
the dynamics. For sampling periods with T different from 
1 sec, the results of Fig. 6 are applicable with the CNRs 
scaled by T. For instance, with T = 10 sec (typical for 
slower processes), the normalized CNR of 0 dB in Fig. 6 
will correspond to the actual CNR of —10 dB. 

It should also be noted that for a given scenario, in- 
creasing A to such an extent that the algorithm cannot 
track the dynamics will lead to degraded system perfor- 
mance. This can be clearly seen from Fig. 6, where it 
is observed that the array gain for A = 0.975 is approx- 
imately 1 dB less than for A = 0.925. Finally, it can be 
observed from Fig. 7 that in contrast to the array gain re- 
sults, there is little difference in rms phase error between 
the least-squares and center-feed outputs. 

B. Simulation Results for Least-Squares Algorithm II 

Figures 8(a) and (b) plot the performance of least- 
squares algorithm II in terms of signal reconstruction. In 
these figures, the dashed graphs depict the amplified signal 
amplitude at the output of various feeds, while the solid- 
lined graphs represent the amplitude of the reconstructed 
field, i.e., the magnitude of various components of h(k ) 
obtained by the circular convolution of the weight vector 
/(Jfc) obtained from Eq. (19) or Eq. (20), and the received 
signal vector g(k) ) for two different time indices equal to 
30 and 185, respectively. As for the case of least-squares 
algorithm I, it is assumed that the weight vector / is nor- 
malized to have its norm equal to 1. This ensures that the 
noise variance at various points of the reconstructed field 
is equal to the variance of the input noise field and, thus, 
the comparison in terms of signal amplitudes is equiva- 
lent to comparing the reconstructed signal-to-noise power 
ratio. The results in Fig. 8 correspond to the same set 


of signal parameters as for least-squaTes algorithm I and 
a CNR of 10 dB. As is apparent from Fig. 8, the least- 
squares algorithm II focuses most of the signal power that 
is originally dispersed in 16 taps into the center tap. A 
more appropriate measure of the effective focusing is ob- 
tained by the signal amplitude of the center tap of the re- 
constructed field, which is plotted versus time in Fig. 9(a). 
In this case, only -5.2 dB of the received signal power 
is scattered in the other taps for the reconstructed field. 
Figure 9(b) plots the phase error of the center tap signal 
and shows an rms phase error of 0.12 rad. Figure 10 plots 
the corresponding results for the case of 7o = 100 and dif- 
fers insignificantly from the corresponding results in Fig. 9. 
Thus, the least-squares algorithm II simultaneously opti- 
mizes the signal-to-noise ratio. Figure 11(a) plots the sig- 
nal loss in the center feed of the reconstructed signal for 
the least-squares algorithm II. As shown in the figure, with 
the parameter A = 0.925, a loss of 1.25 dB can be achieved 
for the high CNR case, which is similar to that obtained 
for modified least-squares algorithm I. Figure 11(b) plots 
the corresponding results for the rms phase error of the 
reconstructed signal, showing an rms phase error of about 
0.12 rad at CNRs higher than 10 dB. 

Figure 12 plots the results when the algorithm is ap- 
plied to a configuration of feeds connected to amplifiers 
with different noise figures. For this example, the case 
wherein one-half of the total number of amplifiers have 
6.0 dB higher noise temperature than others is considered. 
The CNR in the figure is still measured with reference to 
the amplifier with the lower noise temperature. As may be 
inferred from Fig. 12(a), the asymptotic signal loss (CNR 
> 10 dB) in this case is 2.2 dB as opposed to 1.2 dB for the 
case in which all of the amplifiers have low noise temper- 
ature, thus resulting in only 1 dB additional degradation. 
Note that if all the amplifiers were replaced by ones with 
higher noise temperature, then the degradation would be 
about 3 dB with reference to this lower CNR, thus result- 
ing in an effective loss of 9 dB. 

It may be remarked that in the above presentation, the 
sampling period T has been normalized to 1 sec, but the re- 
sults are also applicable to different sampling periods by a 
simple normalization. As is apparent from Eq. (21), while 
increasing the sampling rate by a factor K t one should cor- 
respondingly reduce the actual CNR by the same factor to 
obtain the algorithm’s performance for this case. 

VII. Conclusions 

From the simulations presented in the article, it can 
be observed that for the relatively fast distortion process 
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with a moderate dispersion and the array geometry con- 
sidered, the multi-element array configuration provides an 
improvement of about 3 dB over a single-feed system. For 
a slower process with possibly higher spatial dispersion, 


the improvement is expected to be higher and to a certain 
extent will also be influenced by a match between the ar- 
ray geometry and the pattern of the received signal power 
dispersion. 
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FOCAL PLANE FEEDS 



Fig. 1 . RF signal combining In focal plane. 
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Fig. 2. Least-squares algorithm I: (a) reconstructed signal amplitude; (b) reconstructed 

signal phase. 
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Fig. 5. Least-squares algorithm I: (a) signal combining loss versus 
CNR for reconstructed signal; (b) rms phase error versus CNR for 
reconstructed signal. 
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Fig. 6. Signal combining loss versus CNR for reconstructed signal: (a) A = 0.925; (b) A = 0.95; (c) A = 0.975. Also (d) modified 

least-squares array gain versus CNR. 
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RMS PHASE ERROR 



Fig. 7. rms phase error versus CNR for reconstructed signal: 
(a) A = 0.925; (b) A = 0.95; (c) A = 0.975. 




Fig. 8. Focal plane reconstructed field amplitude, least-squares 
algorithm II: (a) after 30 samples; (b) after 185 samples. 
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Fig. 10. Least-squares algorithm II, 7 0 = 100: (a) focal plane reconstructed signal 
amplitude versus time; (b) focal plane reconstructed signal phase. 
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Fig. 12. Modified least-squares algorithm I, different noise fig- 
ures: (a) signal combining loss versus CNR for reconstructed field; 
(b) rms phase error versus CNR for reconstructed field. 







Appendix 


The following shows that the modified least-squares 
algorithm I of Section IV achieves constrained maximiza- 
tion of the signal-to-noise ratio. In the first instance, the 
time averages are replaced by the ensemble averages. 

Denoting by $,(&) the signal component of the ith 
array element output, consider the problem of minimizing 


E |j5| 2 j but only one of these would simultaneously min- 
imize the error function, Eq. (A-l). It may be remarked 
that there are, in general, an infinite number of solutions 
that satisfy Eq. (A-4), and in effect these orthogonalize the 
estimate S and the “error” (5—5). Among these solutions, 
the one maximizing E |j5| 2 j is selected. Equation (A-4) is 
termed the orthogonality constraint. 


H = E 


N 




Wi Si I 


i= 1 


(A-l) 


with respect to w iy i = 1,2,... N. Setting the partial 
derivative of H w.r.t.Wj to zero yields 


E 


N 


S-J2 


5=1 


(A-2) 


Now with 5 ^ w H s. and s ^ [s \ , 52, . . . , s^Y » index H 
may be written as 

E [|5 - 5| 2 ] = E [|5| 2 + |5| 2 - 55* - 5* S] (A-3) 


At the optimal point, the following is obtained from 
Eq. (A-2). 


E [(5 - 5)5*] = 0 (A-4) 


The optimization index, Eq. (A-l), does not include 
the noise variance at the output of the array combiner, 
which is given by cr 2 ||u;|| 2 where cr 2 is the variance of 
i/j(ib), the noise at the input of the combiner. Thus, now 
Eq. (A-l) is minimized subject to the constraint 

INI 2 = K (A-6) 


for some constant K . Or, one can simply minimize 

E[\S-u,»s\ 2 ]+l3(\\w\\ 2 -K) (A-7) 

where j3 is the Lagrangian multiplier. An analysis similar 
to derivation of Eq. (A-5) shows that with a constraint 
similar to Eq. (A-4), the index is given by 

\S\ 2 - e[\S\ 2 ] - 2/3K (A-8) 


for some constants /? and K . Thus, again the algorithm 
maximizes E ||5| 2 j subject to the constraint that the out- 
put noise variance is equal to a constant Kcr 2 , and thus 
effectively maximizes the output signal-to-noise ratio. 


Adding the left-hand side of Eq. (A-4) and its complex 
conjugate to Eq. (A-3) yields the following form for the 
optimization index, subject to the constraint, Eq. (A-4). 

H = |5| 2 - E [|5| 2 ] (A-5) 


Thus, the algorithm that minimizes Eq. (A-l) also maxi- 
mizes E ||5| 2 |, subject to the constraint, Eq. (A-4), i.e., 
it is also a signal maximization algorithm. There may 
also exist solutions that optimize E |5| 2 without the con- 
straint, Eq. (A-4), but these may result in large phase error 
with respect to 5, the desired signal, i.e., 5 and Se** (for 
any random phase <j>) both have the same value of the index 


Now from the independence of the received signal s* 
and noise it follows that 

E [|5 - w H i\ 2 } = E [|5 - w H s | 2 ] + HI 2 * 2 (A-9) 

Minimization of Eq. (A-9) subject to the constraint, 
Eq. (A-6), is thus identical to the minimization of Eq. (A- 
1) subject to the constraint, Eq. (A-6), and thus effec- 
tively maximizes the signal-to-noise ratio under the near 
orthogonality constraint. Now for the large value of k, 
the index k~ l Jk with J* given by Eq. (3) approaches the 
left-hand side of Eq. (A-9) under appropriate ergodicity as- 
sumptions, and the algorithm of Section IV thus achieves 
constrained optimization of the signal-to-noise power ratio. 
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It may be remarked that the least-squares algorithm 
in the absence of the constraint, Eq. (A-6), effectively min- 
imizes the following index 


|Sf-£ |Sf 


<t 2 \\w \\ 2 


(A-10) 


subject to a constraint similar to Eq. (A-4). Since the last 
term in Eq. (A-10) represents the noise power at the com- 
biner output, it can be observed that the standard least- 
squares algorithm effectively maximizes the sum of the sig- 
nal plus noise power rather than the signal-to-noise power 
ratio. 
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